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We discuss the conclusions on the symmetry of hidden order (HO) in URu2Si2 that may be drawn 
from recent torque experiments in rotating magnetic field by Okazaki et al.— . They are very sensitive 
to changes in the magnetic susceptibility induced by HO. We show that the observed twofold angular 
torque oscillations give evidence that hidden order has degenerate E- type (yz,zx) symmetry where 
both components are realised. The oscillations have the wrong characteristics or are absent for the 
ID nontrivial representations like quadrupolar Bi{x^ — y^) and B2{xy) type HO or hexadecapolar 
A2{xy{x^ — j/^)) type HO. Therefore they may be excluded as candidates for hidden order. We also 
predict the field-angular variation of possible field-induced Bragg peaks based on underlying E-type 
order parameter and discuss the expected elastic constant anomalies. 
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I. INTRODUCTION 

The moderate heavy fermion compound 
URu2Si2 serves as a paradigm of electronic hid- 
den order (HO) and field induced quantum phase 
transitions. Although investigated by many techniques 
since considerable time its microscopic nature and in 
particular its symmetry remains elusive. Progress in 
understanding was greatly advanced through pressure 
experiments which demonstrated phase separation of 
hidden order (HO) phase and parasitic heterogeneous 
antiferromagnetism (AF) with tiny moments (see Ref. [3 
and references cited therein). At low pressure the 
bulk of high quality samples exhibits only HO with 
ordering temperature T/i=17 K which coexists with 
the superconducting phase appearing below Tc =1.3 
K. The HO transition at Th is associated with a large 
entropy release that can be accounted for by a dramatic 
increase in the gap of incommensurate itinerant spin 
excitations^. In addition to ARPES'' and STM^ and 
transport''"^ experiments this shows that a microscopic 
model for the hidden order has to take into account the 
partly itinerant character of 5f electrons which has been 
discussed extensively in Ref. The microscopic models 
considered sofar are too numerous to be discussed here 
(see references cited in Ref. l). The purpose of this 
work is a more modest one: To reconsider the question 
of proper symmetry class of the HO parameter from 
analyzing new oscillatory torque measurements^ and 
previous elastic constant measurementai^. We will 
propose that the former lead to a unique answer for HO 
symmetry in URu2Si2. 

For the analysis of symmetry classes it is conve- 
nient to start from localised 5f electron states which 
can carry multipolar moments up to rank five in 5f 



angular momentum operators. This approach allows 
a systematic and physically appealing classification of 
HO symmetry. Each of those with rank > 2 may be 
candidate for hidden order because unlike the dipolar 
(rank 1) moments they are hidden in zero field neutron 
diffraction experiments without magnetic field. The 
full classification is given in Table H Indeed the first 
serious discussion of HO symmetry in URu2Si2^^ used 
such localised picture. The tetragonal symmetry group 
D^fi of URu2Si2 has three nontrivial ID representa- 
tions (A2, Bi and B2) and one doubly degenerate E 
representation. Any possible multipolar HO has to 
belong to one of them denoted by F. Since pressure 
dependence shows that HO is in competition with AF 
order at the commensurate Q — (0,0, 1) wave vector (in 
bet structure this is equivalent to Q = (1, 0, 0)) one may 
expect the same (antiferro) ordering wave vector for HO. 
This is supported by the observation that the resonance 
excitation found at Q by inelastic neutron scattering 
(INS) turns into the elastic AF peak above the critical 
pressure of HOi2. and by the similarity of Fermi surfaces 
in the AF and HO phases^'^. Furthermore the absence of 
any macroscopic lattice distortion^** below is evidence 
against the ferro HO case with Q = 0. Nevertheless we 
will also include ferro-type HO in our analysis. Most 
frequently quadrupolar HO type (Or) has been invoked. 
It was proposed to be of Bi(02,2_j^2) typei^ or B2{0xy) 
typeiiii^ii^. In addition higher order multipole moments 
like octupoles of B2{Txyz) typeiiiii and hexadecapolar 
order A2{llxy{x^-y2) y^^^^ have been suggested and 
discussed. For a general review of multipolar type order 
in crystalline electric field (CEF) split f-electron systems 
see Refs. [HM 



The most powerful means of identifying the HO con- 
sists in inducing dipolar (magnetic) order at the same 
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wave vector by application of an external field. The in- 
duced magnetic order can then be seen in neutron diffrac- 
tion. The directions of induced moment as function of 
field orientation reflect the symmetry of the underlying 
hidden order which can then be identified. For example 
by this method the HO in some rare earth skutterudite 
compounds has been determined^Si^i. The HO should 
also lead to subtle changes in the uniform dipolar (mag- 
netic) and quadrupolar susceptibilities. The latter can 
be probed with high sensitivity through the elastic con- 
stant anomalies below the HO transition. The magnetic 
susceptibility tensor in the HO phase has recently been 
investigated with the equally sensitive torque method^ 
below Th which we will discuss first. Above the HO 
transition the susceptibility tensor should have fourfold 
symmetry in the tetragonal ab plane corresponding to 
fourfould oscillations in the torque. Below Th this may 
no longer be true and twofold torque oscillations may ap- 
pear depending on the symmetry of the HO. Their obser- 
vation therefore allows to draw conclusions on the latter. 
The torque exerted on the sample of volume V by an ex- 
ternal field is given by t = iiqVM. x H with M — xH 
where x is the susceptibility tensor. For field rotation 
in the ab-plane with H = H {cos if, simp) the non-zero 
component Tz{(p) = T{ip) is given by 



^ {Xxx - Xyy ) sin 2^3 - Xxy COS 2ip 



(1) 



where (p is counted from the a axis and Xxy — Xyx is 
assumed here. To analyse the torque oscillations it is 
therefore necessary to compute the susceptibility tensor 
above and below the HO transition. 

In Sec. iniwe first give a symmetry classification of pos- 
sible multipole order parameters in URu2Si2 within a lo- 
cal f-f coupling scheme, in analogy to Ref. [221. In Sec. IIIII 
a Landau theory for selected even rank multipoles is used 
to calculate the torque oscillations. In Sec. IIVI we con- 
sider the possibility of field induced staggered moments 
in the HO phase and in Sec. |V] we discuss the anomalies 
of some symmetry elastic constants at the HO transition. 
Finally Sec. IVII gives the conclusion. 



II. SYMMETRY CLASSIFICATION OF 
MULTIPOLAR ORDER PARAMETERS FOR 
f-ELECTRON SYSTEMS 

In /-electron systems, due to strong spin-orbit cou- 
pling, /-electron states are described well by the lower 
j—5/2 multiplet where j is the total angular momentum. 
The multiplet consists of three Kramers doublets. Then, 
due to six states, irreducible tensors for representing mul- 
tipoles are available up to rank 5. With increasing rank 
k they are denoted as dipole (J, rank 1), quadrupole (O, 
rank 2 ), octupole (T, rank 3) hexadecapole (H, rank 4), 
and dotriacontapole (D, rank 5) moments. The latter 
will not be included in the discussion and the tables. 



The irreducible tensor has 2fc-|-l components 
Jq'^'^ < which satisfy the following relations 



q 



(2) 
(3) 



where J±= ± ijy are raising and lowering operators, 
J" is a-component of total angular momentum operator. 
When these multipole moments are in a CEF produced 
by cubic or tetragonal surroundings, the corresponding 
tensors have to be classified according to irreducible rep- 
resentations of the point group^ which is here. It 
is convenient to express the multipole operators that act 
on CEF states by single particle operators adapted to the 
-D4/1 symmetry. We start from the tetragonal CEF states 
of U in f-f coupling scheme given by 



±) = \7 = l,±)=a\±-)- yr^l T ^ 



IF^') , ±) = |7 = 2, ±> = v/T^I ± ^) a| T 2 ) 
|r6,±) = |7 = 3,±) = |±i) 



(4) 



where a is a real number determined by the parameter set 
of tetragonal CEF. Furthermore, a denotes pseudo-spin 
states of /-electrons in the Kramers doublets which are 
connected by time reversal and 7 the various doublets. 
Two Kramers doublets of 7=1 and 7=2 correspond to Fy 
irreducible representation in the tetragonal crystal struc- 
ture while the remaining one with 7=8 belongs to Fg 
irreducible representation. 

By using these CEF bases for /-electron states, we 
rewrite operators of irreducible tensors Jq'^K When q is 
an even (odd) number, the matrix element (7(t| Jg'^^ |7'ct') 
can only be non-zero for a— a' {(j=—a') Using this prop- 
erty irreducible tensors with q=2p are given by 



(5) 



7,7' 



where the matrix element (70'|</2p^ I7'''') satisfies 

(7^l4,"l7V) = (-)'=(7'^l4p^l7^) 

= {l'<^\J%M (6) 
Similarly, irreducible tensors with q=2p + 1 are given by 



4';ii - EE(^'^i'^2plii7''?)4./7'. 



with 



(7a|jgiil7''?) = (-)'+'(7V|jglil7'?) 



(7) 



iy^\J%,+i)M (8) 
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TABLE I: The first column shows the irreducible representa- 
tion of D4h. The second and third ones describe corresponding 
bases consisting of multipole moments and their one-particle 
operators (Eq. (O), respectively. Here, in the second column 
of A^ irreducible representation, Nf is the number operator 
of /-electron. The superscript ± of irreducible representation 
expresses the parity with respect to time reversal. 



of Table HI The hidden order parameter must belong to 
one of these representations. In particular they may be 
used to construct microscopic models for HO. Although 
this is not the topic of the present investigation we briefly 
describe the principal type of possible microscopic mod- 
els. 

In primarily localised models the multipoles of CEF 
spli^iii 5f states given in Table H] are coupled by effec- 
tive intersite interactions obtained from eliminating the 
conduction electrons, see e.g. Ref. I2I. Depending on the 
leading interaction channel and appropriate symmetry of 
lowest CEF states a specific multipole order belonging 
to one of the representations in Table U may be selected 
as the ground state. Although this is a local picture it 
should be noted that staggered order of multipoles may 
have strong feedback on the itinerant states by gapping 
large parts of the Fermi surface. This mechanism has 
been invoked in Pr- skutterudites and has been proposed 
as similar in nature to the observed carrier depletion in 
the HO phase of URu2Si220. 

In the itinerant models the hidden order parameter re- 
sults from a pairing mechanism of electrons and holes 
leading to an unconventional density wave or 'nematic' 
state. In such models the gap function (either in spin 
singlet or triplet channel) breaks the crystal symmetry 
and has a nontrivial (e.g. quadrupolar) type form fac- 
tor or k-dependence. These models have been frequently 
proposed as candidates for hidden order in itinerant sys- 
tems, e.g. in Refs. i25^2& . They naturally lead to a partial 
gapping of the Fermi surface in the HO state. Such exotic 
pairing states can only be stabilized in extended Hubbard 
type models for the itinerant 5f bands. 

In this work we do not proceed to the construction 
of a microscopic model for HO but rather try to draw 
conclusions from a symmetry analysis of predicted torque 
oscillations based on Landau theory of HO. 



III. LANDAU THEORY OF TORQUE 
OSCILLATIONS IN THE HO PHASE 



Using these expressions, all Jq finally are given by 
simple linear combinations of charge operator p^ji and 
pseudo-spin operators .S*"^, (a=x, y, and z) defined by 

cr a, a' 

where f " are the Pauli matrices and /^g. creates f elec- 
trons in CEF Kramers doublet state I7, cr). These charge 
and pseudo-spin operators which are either symmetric or 
antisymmetric in the orbital space are convenient as new 
bases of multipole moments. Thus, we introduce a new 
Hermitean basis {<fi^} which is shown in the third column 



As candidates for HO we consider only the even repre- 
sentations of Table |T] because rank 1 dipoles as primary 
order parameter are obviously excluded and rank 3 oc- 
tupoles may induce dipoles as secondary order even at 
zero fielctii. In particular we discuss rank 2 nondegener- 

B+ B+ 

ate quadrupole order parameters (p = (pi'- , (p — (pi^ or 

twofold degenerate {(j)a,(f>b) — {(Pix T'Piy )• I^i addition 
we consider the rank 4 nondegenerate hexadecapolar or- 

A+ 

der (f> = (j>^^ . The Landau free energy functional F con- 
tains the uniform magnetisation M and the HO param- 
eter. For the ID representations it can be constructed 
from invariants of the D^h group using only terms up to 
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fourth order in the HO and uniform magnetization: 

F%o{4>} = act>^ + \hcl^^ 

Bi : Fnom = Fho{<P} + 

gml - Ml) + ~gct,\Ml + M^) 

B2 : Fho{4>} - F%o{4>} + 

g4)M,My + .#2 + M^) (10) 

A2 ■ Fho{<I^} = Flo{<t>} + 

gct^M^MyiMl - Ml) + gcf^^Ml + M^) 

Note that the terms ~ g which are hnear in the order 
parameter (p appear only for ferro- type HO and are 
strictly forbidden (172 = 0) for antiferro-type HO, e.g. 
Q = (0,0,1) due to the requirement of translational 
invariance. In this case only the terms ~ (j)^ contribute. 
Because for all ID HO the square of (j) is an invariant 
it can only couple to the modulus of the magnetisation 
which does not break the in-plane symmetry. The 
argument would also apply for even terms (the only ones 
allowed for Q ^ 0) of higher than fourth order. This 
fact has important consequences for the possible torque 
oscillations in the ID HO phases. 

For the twofold degenerate E representation the free 
energy is 

F°Ho{^aAb} = a{4>l+cf>l) + h,{4'a+<Pt)+b2fyl 
FH0{4>a.4>b} = F''Ho{^a,M+ga{4>l+^l){Ml + Ml) + 

g^[<g - cj>l){Ml - Ml) + 2g2cl)a4>bM,My 

(11) 

In Eqs. ([TUITTt we have a(T) = ao{T - T,-,) and 
fl0i^j^i,2 > 0. Finally the purely magnetic part of the 
free energy given by 

F„,{M} = — + -Xi{M^ + M^) + X2M^Ml^M-U 

(12) 

has to be added. Here H= {Hx,Hy) and M = {Mx,My) 
are external field and uniform magnetisation respectively 
lying in the tetragonal ab plane. The magnetization is 
determined by M = xH. In zero field the susceptibility 
tensor is diagonal with Xxx = Xyy — Xo equal to the 
background susceptibility xo of URu2Si2 which is almost 



temperature independent for field in the plane--. In an 
external field in general Ax = Xxx - Xyy and Xxy,Xyx 
may become different from zero due the underlying 
hidden order. 

At zero field minimization of the free energies for ID 
representations leads to the order parameter {T < Th) 

02^^(T,-T); Fo = -i|(T,-r)2 (13) 
For degenerate E representation there are two possible 
solutions characterised by the vector cj) = [(j)a,4>b) = 
(j){ria,'i]h). Because </> is a commensurate (Q = or 
(0, 0, 1)) order parameter of multipolar density type the 
amplitudes {4>a, 4>b) may be chosen as real. The two possi- 
ble phases E{ria, rjb) are then determined by the numbers 
ria = 0, ±1 and rib = 0, ±1 according to 

E{0,1) ih>b2): 

02 = ^(T, _ T); Fo = -l^in - T)' (14) 
E{1,1) {bi<b2): 

Oi +02 Oi + 62 

All other combinations of E{T]a,rii,) represent different 
domains of these phases. There are four domains with 
equal Fq in each case: E{0,±1), i?(±l,0) for hi > 62 
and £;(±1,±1) for hi < 62. 

For finite field the free energy functional has to 
be minimized simultaneously with respect to the HO 
parameters and the induced uniform magnetization. 
After some lengthy but straightforward algebra the 
solution for all ID and 2D representations may be 
obtained. From this we calculate Ax = Xxx — Xyy and 
Xxy^Xyx as function of temperature and field for each 
symmetry which then determines the torque oscillations 
according to Eqs. (I1I15|) . We give the set of results 
for all HO symmetries of interest and then compare to 
experimental observations to draw conclusions on the 
prefered symmetry of HO. 

The total torque in Eq. ^ may always be split into 
contributions with twofold and fourfold oscillations de- 
noted by t!^"^ and Ti^"*-* respectively. The total T{ip) may 
be written in reduced form as 



r(2) 



r(4) 



T(ip) = To{H) —{Xxy 



Ax 1 

Xyx) + 7r- sin 2<^ - {Xxy 

I 



Xyx)cos2ip 



(15) 



with the scale of the torque given by tq{H) = iiqVH^xo- 



In this slightly more general form of Eq. ([T]) it is taken 
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into account that HO may possibly break reflection 
symmetry with respect to [110] diagonals in the plane, 
leading to Xxy 7^ Xyx for general field direction. Below 
we compile the results for the torque oscillations for the 
main candidates of HO symmetry. This is most conve- 
niently done in a reference frame where the field direction 
refers to the crystal axis, i.e. H = H{cos(p, sin(p) with ip 
counted from a-axis or [100] direction. Experimentally 
however the oscillations are preferably counted from the 
[110] direction i.e. with H = ff(cos(<^ + f ),sin((^ + |)) 
with If = + J. For clarity we will give two- and four- 
fold oscillations for both conventions. They are simply 
related by using the identities sin2(^ = cos2ip, cos 2^3 = 
— sin2(^, siiiAip — — sin4(^, cos 4^3 — — cos4(p. For 
comparison with the experimentally observed oscillation 
pattern we will use the (p convention . 

First we discuss the results for the nondegenerate 
HO phases. For convenience we introduce the constants 
5o = 2xl{X2 - Ai) and Aq = gxoiaofh/b)-2 . The fully 
symmetric terms in the Landau functional Eq. (jlOp 
only lead to a shift of the HO transition tempera- 
ture in the field according to Th = Th — AT/^ with 
AT/i = (gXo/ao)(xo-ff^)- Then we obtain for 

Bi{x^ — y'^) quadrupolar HO: 



B2{xy) quadrupolar HO: 



Ax 

Xo 

Xxy 

Xo 

V 
^(4) 



SoixoH^)cos2ip 

ToAo(l —)2 cos2(^ 

Th 

■^ToS ixo H'^) sin Aip 
-roAo(l - sm2ip 
-^To^o(xo^^)sin4(^ 



(17) 



This model also has both twofold and fourfold oscil- 
lations. Again for HO with Q ^ the former vanish 
identically due to g, Aq = 0. Furthermore for the 
case Q = where they are non-zero the amplitude 
increases ~ (T — Th)^ as in the Bi case. This is a 
much too strong increase below Th in comparison with 
experiment that seem to follow a weaker dependence 
Tf{H,T) ~ ci(T - fh) + C2(T - fhf. Therefore the 
B2{xy) HO model cannot explain the observed torque 
oscillations in URu2Si2. 

A2{xy{x^ — y^)) hexadecapolar HO: 



Ax 
Xo 

Xxy 
XO 



r(2) 



-4A, 



o[l-^)' +So{xoH^)cos2ip 



(16) 



r(4) - 



T 1 

-2roAo(l- ;^)3sin2^ 

^T-o<5o(xo-ff^)sin4i^ 
T 1 

— 2roAo(l cos2ip 

Th 

--^ToSQixoH"^) sin Aip 



This model in principle leads to both twofold (period 
7r) and fourfold (period tt/2) torque oscillations with 



amplitudes 



{H, T) and 



respectively. 



However for lio with Q 7^ as is presumably the case 
we have 5, Ag = and therefore rj^^ (H, T) = , i.e. no 
twofold oscillations appear. Furthermore even for Q = 



(2) 

when t\ '{H,T) is nonzero their angular dependence 
~ cos 2ip has the wrong phase with maxima and minima 
interchanged compared to experimetal observation. 
Therefore Bi{x^ — y'^) type quadrupolar HO is not 
compatible with torque experiments. 



1 

2x^ 
1 

2x^ 



{Xxy 
iXxy 



Ax 

Xo 

Xyx)) 

Xyx)) 

V 

V 



So{xoH'^)cos2ip 

-2(xoAo)(xoi?') 

-(xoAo)(xoH')( 



-^ro{xoH'^)[So sin4(/3 



T 

't\ 

T 



cos 2ip 



(18) 



+4(xoAo)(l - ^ 

J^h 



COS iip] 



-^7-o(xo-ff^)['^osin4^ 
-4(xoAo)(l-;~p)'cos4(^] 

J^h'^ 



The A2 symmetry is not realized for quadrupolar 

(rank 2) HO because it only breaks reflection symmetry 

with respect to [110] type planes but not C4 rotational 

symmetry. It appears first as rank 4 (hexadecapolar) HO 

(2) 

parameter. However in this case = 0, independent 
of Q , in conflict with experiment. Furthermore for 
Q = (Ap 7^ 0) the fourfold oscillations below Th would 
be skewed with respect to the [110] diagonals due to 
the abovementioned lack of reflection symmetry of the 
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hexadecapolar HO. This feature is also incompatible 
with torque experiraents. We conclude that the hexade- 
capolar order parameter cannot explain the observed 
twofold torque oscillations. 

In summary if we restrict to considering only Q 7^ 
HO then none of the ID representations can produce the 
twofold torque or susceptibility oscillations that have 
been observed by Okazaki et al^ below Th- For Q = Bi 
and B2 quadrupolar HO may lead to twofold oscillations 
but they have the wrong temperature behaviour or 
angular dependence and therefore are not compatible 
with experiments. Furthermore the A2 hexadecapolar 
OP will never lead to the observed twofold torque 
oscillations. The failure of all ID HO parameter models 
cannot be remedied by including higher order terms 
(in 0, M) in the Landau functionals. They would only 
renormalize the amplitudes but would not change their 
principal characteristics. 

Therefore we are lead to consider the twofold degener- 
ate E representation with its two different phases as only 
remaining candidates. In this case the fully symmetric 
terms in Eq. (|lip lead to a shifted transition temperature 

fh = Th- ATh with ATh = {goXo/ao){xoH'^)- Then we 
obtain for 

E{0, 1)(0, zx) quadrupolar HO: 



Ax 

Xo 

Xxy 

Xo 
Ao 

^(2) 

rW 

V 



.(4) _ 



Th 







(<7iXo)(aoT,/6i); 5o = 2xg [(A2 - Ai) + ^ 

L Oi 

T 

2toAo(1 - -^)sin2v9 (19) 
Th 

■^toSq{xoH^) sin Aip 
T 

2roAo(l - -^)cos2.^ 
Th 

-^To'5o(xo-ff^)sin4(^ 



The expressions for the E{1,Q) domain may be 

obtained by substituting Aq — ^ — Aq in the above 

expressions. It means that twofold amplitude changes 
(2) 

— J> —T^ when we change from E(0,1) to E(1,0) 



(2) 

Sign t\ 
domain. 



This is the first model which leads to twofold os- 
cillations r^^ 7^ even in the case Q 7^ 0. This is 
because gi may be finite also for finite Q since the 



quadratic term ^ (0^ — 
invariant. However t1 



in Eq. ([TT]) is translationally 
~ cos2(^ and therefore its 



correct experimental positions. Due to these wrong 
angular characteristics the E(0,1) or E(1,0) phase can 
also be excluded as the HO candidate. Finally we discuss 

^^(l, V){yz, zx) quadrupolar HO: 



— = 5o(xoi^')cos2v? 



Ax 

Xo 

Xxy 

Xo 
Ao 

^(2) 

V 



-2A, 



92Xo{aoTh/ibi + b2)) 



2x 



(A2 - Ai) + 
T 



2toAo(1 — ^)cos2iy9 
Th 

^T-o^o(xo-ff^)sin4(/7 
T 

-2roAo(l - -^)sin2(^ 
Th 

-7'^o^o(xo-ff^)sin4(^ 



(20) 



This model has the proper behaviour of torque angu- 



lar dependence r^' 



(2) 



, T^''^ in comparison to experiment as 
shown in Fig. [T] The fourfold oscillation is present al- 
ready above Th and temperature independent in the lead- 
ing order. On the other hand the twofold oscillation sets 
in below Th and grows linearly with (l — T/Th) in leading 
order. If we also include terms ^ gi{4>t ~ 4'b)i^x ~ ^y) 



the temperature dependence will be given by 



Ao(l-^) 



Ao(l — 7^)^. It appears that in experiment 

^ h 

both terms are present, at least in the samples with small 
size. Experimentally it was also found that rj^'' , rj^^ have 

opposite signs which implies A2 < Ai and So < 0. We 
conclude that the basic features of the torque experi- 
ment in Ref. — are reproduced by the E(l,l) quadrupolar 
HO model. This is in agreement with the order param- 
eter proposed from analysis of resistivity anisotropy in 
Ref. [2^. As explained above there are four domains with 
iVajVb) = (iljil) in this phase. The domain type mat- 
ters only in the last term of Eq. (ITT|) through the sign 
oi rjaTib ~ ±1- For domains with r\a"']h — ~1 one has to 
replace Aq — > — Aq in Eq. ((20)) leading to opposite sign 
for the twofold oscillations. If domains with ^a^h ~ ±1 
are equally distributed their twofold amplitudes will can- 
cel. Indeed the latter are most prominent in the smallest 
samples where one would expect a single domain phase 
like E{1,1) to be realised. Finally we may rewrite the 
twofold oscillation in Eq. (pn|) in the form 



A'' 



(2) 



maxima/minima are shifted by 7r/4 with respect to the 



Af sin 2^ 

l/(^)Mo = ^[M[no] 
Xo 2 ' ^ 



Mriiol] (21) 
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where Mq = XqH- This shows explicitly that the 
twofold amplitude in the E(l,l) phase is due to the 
off-diagonal susceptibility induced by the background 
HO. The last identity may be obtained from M{ip) ~ 
{XoH)[l + {xxy/xo) sm2ip]. It signifies directly the break- 
ing of fourfold symmetry in the HO phase with respect 
to the diagonal [100] and [110] axes. The linear relation 
in Eq. (|21|) is found experimentally for low fields. For 
larger fields there is however a sign reversal. Since the 
latter depends on the type of domain we speculate that 
this may be due to change of domain preference at higher 
fields. 



r(C4v) O^^ (iC iEYen) O^' 0^ (Odd) 

Fa Ol <l>f' = ;^(P23 + P32) Tj ^ ^{Sh + Sh) 

T4 Oxy 4>1^ = TfC-gfa - ■S'm) r^HZ <^f' =7g(P23-p32) 

I5 Uyz (p\x ~ ^\Ol2 ~ J2l) Jx Cplx ~ Jll 

Tra J^E'^ i fox ox \ rj^a ox 

J^x V'ix — "yf 1'523 ~ '-'32j J^x 92x —'322 

Hx 4>ix = "Tg^*^^^ ~ ■S'13) ~ ■S'33 

Ozx (t>fy ~ ^{3^2 ~ S21) Jy 4'ly = 'S'll 

Hy 4>2y = ■^(S'la — 5*32) Ty (f2y ~ S22 

(fliy = -"Tjl'S'll — •S'13) (t)§y — S'33 



IV. CLASSIFICATION OF POSSIBLE FIELD 
INDUCED DIPOLE MOMENTS 

We have shown that the results of torque experiments 
and their Landau theory analysis strongly favor the an- 
tiferroquadrupolar HO with E(xz,yz) symmetry of the 
(1,1) type and wave vector Q = (0,0,1). One of the 
most direct ways to look for hidden order is to identify 
the direction of induced staggered dipolar moments in an 
external field. For field along [001] or [100] directions the 
symmetry is reduced from D4/1 to 04^, or C2V respectively. 
The proposed candidates for primary underlying hidden 
order (quadrupolar or hexadecapolar) have all even time 
reversal parity. Then a magnetic field may induce dipo- 
lar or octupolar moments which have odd time reversal 



TABLE II: The first column shows irreducible representations 
of C4v for magnetic field along [001] direction. The second and 
third ones describe corresponding bases consisting of multi- 
pole moments and one-particle operators with even parity, 
respectively. Those with odd parity are given in the fourth 
and fifth columns. 



parity. The full compatibility tables for the induced mo- 
ments and field along the two symmetry directions are 
given in Tables |ll] and IIIII Note that more cases of HO 
symmetry have been included in the tables as compared 
to the previous Landau analysis. Table ITTl shows that for 
field along [GDI] the E(l,l) type quadrupolar HO may 
have induced moments along the [110] diagonal direc- 
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tions. It is also obvious from Table lllll that for quadrupo- 
lar E (0,1) (O^^) type HO and field along [100] the in- 
duced staggered moments mq should be oriented along 
[001] direction. The same holds true for E(1,0) (Oyz) 
and field along [010]. If the field is perpendicular to the 
domain orientation of HO no dipolar moment but only 
Txyz octupolar moments may be induced. For the E(l,l) 
phase both quadrupolar components lead to an induced 
moment along [001]. This is indeed the same direction as 
observed in the genuine AF phase of URu2Si2 above the 
critical pressure of 0.7 GPa^. The moment dependence 
for general field direction in the ab-plane may be calcu- 
lated from an extended Landau functional including the 
staggered magnetisation {uizq) terms. It is given by 

F - FHo{<l)a,M + F,n{M) + -^ml, 

+>^Q{HxmzQ(f)a + HyUizQcj)!,) (22) 

Here xq denotes the staggered susceptibility and Aq a di- 
mensionless coupling constant. Mininimization of F leads 
to the expression for to^q. The physical important quan- 
tity is the intensity of field induced Bragg peaks which is 
~ [wzqI^. It is given by the following domain dependent 
expressions 



r(C2v) ct>i (Even) (Odd) 

-/"f ^ = Pii 7^ <^f." = sf^ 

O2 4>2^^ = P22 (t>2x = S22 

of <t>-j. ^ = P33 Tf (Pix = S3S 

^2 Oyz 4>lx^ = 71 ("^l^ - 'S'll) Txyz 0i ^ = 7^(^12 - p2l) 

Ts Oxy 01^ = -^(S'l2 — S'll) Jy 4>ly ~ Sfi 

H" 4>l^ ~ ^{Sis — S32) Ty (j>2y =5*22 

Hy 4'2y =-^{Slz~Sl2) <^2^ ^822 

Hy 4>3y = -75 ("S'll — •S'i'3) rf (^3^ ~ 



TABLE III: The first column shows irreducible representa- 
tions of C2v for magnetic field along [100] direction. The sec- 
ond and third ones describe corresponding bases consisting 
of multipole moments and one-particle operators with even 
parity, respectively. Those with odd parity are given in the 
fourth and fifth columns. 



V. ELASTIC CONSTANT ANOMALIES BELOW 
THE HO TRANSITON 



E{l,±l): |m,Q|2 = (AQXQ)'i/'(l±cos2(^),^2 
£;(1,0): |m,Q|2 = (AQXQ)2i?'i(l-sin2(^)</.2 

i^(0,l): |m,Q|2 = (AQXQ)'i/'^(l + sin2(^)</.2 

(23) 

where 02 ^ [a^ / {hx+h2)){fh-T) and ^ ^ (ao/6i)(T',,- 
T) . As most promising candidate for HO the E(l,l) 
phase has been identified. In a diffraction experiment 
the intensities at Q and equivalent points should vary 
according to the first of the above expressions. As for 
the torque, twofold oscillations as function of in-plane 
field angle in the intensity are predicted for a single do- 
main. If domains with riarjb = ±1 are equally popu- 
lated only the average intensity proportional to [jtIzqI^ = 
{XqXq)^F['^4'^ should be observed which does not de- 
pend on the field angle. As in the case of torque ex- 
periments, twofold angular variation of Bragg intensities 
should therefore only be observed in small single domain 
samples. Here xq can be expected to be of order l/Th- 
Then {H/Thj^ sets the scale for the field induced mo- 
ment. Since T/i = 18 K is quite large one may conjecture 
that magnetic moments induced out of the HO phase ap- 
pear only at very large fields. This means the intensity of 
magnetic Bragg peaks at accessible fields is likely to be 
dominated by those of the tiny volume fraction (^ 1% 
according to Ref. [29l ) of the parasitic AF phase which 
are already present at zero field—. 



As mentioned in the introduction the absence of any 
homogeneous lattice distortions below T/r^ suggest HO 
at nonzero wave vector. On the other hand small but 
distinct anomalies in the elastic constants below the 
HO transition of URu2Si2 have been found^°. In the 
same manner as torque experiments probe the magnetic 
(dipolar) susceptibilities the elastic constant anomalies 
probe the quadrupolar susceptibilities associated with 
background HO below Th- As in the torque experiments 
for X in Ref. [H the ultrasonic experiments for the 
symmetry elastic constants are sensitive for the 

changes caused by the onset of hidden order. The 
results obtained in Refj^ have, however, sofar not 
been analyzed in their relevance for the HO symmetry 
questions, only the much larger anomalies due to local 
single ion quadrupolar CEF excitations have been 
disussed. Generally it was found that elastic constant 
anomalies Acr below Th are quite small of the order 
of Acr/cr — 10~^ due to a slight change of slope in 
cr (r) . It is not clear at present whether this change is of 
significance and contains information on the symmetry 
of HO. Therefore we investigate this problem in the 
present phenomenological context. 

To perform such analysis we extend the previous Lan- 
dau treatment for zero field and include the coupling 
terms to the homogeneous elastic strains. They can be 
classified (e.g., in Reffiol) into the 1)4/1 symmetry strains 
er associated with elastic constants cr leading to the elas- 
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TABLE IV: Step-like elastic constant anomalies due to linear strain coupling F'^I'Lho ^ type HO phase for transverse 

Cii — Ci2 and 2c66 symmetry. Note that the anomaly is absent for the prefered HO symmetry E(l,l) which is well compatible 
with experiments^". Here Ac-y and Acs are independent of domain type in both phases. 



phase Ac^(cii-ci2) Aca(2c66) OP 



£;(0,1),S(1,0) -£x(i + 2f^<^2)2<0 2(i;f^>0 = f-; {b^ > b^) 



tic energy 

Fei = \Y.44 (24) 
r 

Here Cp(r) is the elastic background constant in quasi- 
harmonic approximation (containing the anharmonic 
contribution and the effect of localised CEF excitations). 
The total elastic constant cr(T) = c^{T) + Acr(T) 
contains the small effect of the coupling to HO below 
T;i and this part may in principle give additional 
information on the symmetry of the latter. 

For D/^h symmetry there are two one dimensional Ai 
type fully symmetric volume strains and in-plane tetrag- 
onal strains of Bi- type = {l/y/2){exx — £yy) with 
C7 = cii - C12 and B2- type eg = \/2ej;y with eg = 2cee. 
In addition the two dimensional out-of plane strain 
of E-type is given by = V2{ezx,eyz) with Ce=2c44. 
The symmetry strains couple to the HO parameter 
via i) first order (linear in er) strain coupling leading 
to corrections Acr in second order of the coupling 
constant and ii) second order (quadratic in ep) leading 
to corrections Acp in first order of the coupling constants. 

We first focus on the linear strain coupling. For HO 
with Q ^ translational invariance is broken and only 
the square (or generally even powers) of the nontrivial 
ID order parameter can couple linearly to the ID ho- 
mogeneous strains. But for any ID representation F we 
have F X F = Ai. Therefore the linear strain coupling of 
the order parameter can exist only for the volume strain 
. This means that the nontrivial elastic constants (aside 
from the bulk modulus) should not exhibit an anomaly 
at Th for Bi, B2 and A2 type HO resulting from linear 
strain coupling. Then the only case that remains to be 
discussed are the two possible E phases. The linear strain 
order-parameter free energy for this case is 

P^t^HO = 97<^ii4'l - (t>b) + 9i<^i{4>i - (t>b) + gsf-s4>a<t>b 
+gses{4>l + ^l)Mb (25) 



Here we did not include the out-of plane E-type C44 strain 
Ee because it has no linear coupling when Q 7^ 0. The 
anomaly of the elastic constants below Th may now be 
obtained from 

Acrm = E*k.(^); 4. = ^ M 

where F = Ff)+Fei+F[l}_jjQ and F = 7, 5 denote cn— C12 
and C66 elastic constants, respectively. After determin- 
ing fy from the equilibrium equations the calculation 
of Acr is straightforward. The final results are shown in 
Tabic IIVI for the two possible E phases. Generally step- 
like Acr(T)^ 5r anomalies at Th are predicted. However 
as mentioned before experimentally rather slope changes 
in criT) are observed. This may be partly due to broad- 
ening of the transition as in the case of the specific heat 
anomaly at Th- On the other hand slope changes are 
naturally obtained as a result of the quadratic strain cou- 
pling which we will discuss now. For the E-type HO this 
mechanism is described by 

F^^sIho = + + 9'Ml + <\>l) (27) 

where g'-^ are the second order strain coupling constants. 
Contrary to first order there is also a possible second 
order coupling of strains (C44 elastic constant) but it 
will not be considered here. Because ep (F = 7,(5) trans- 
forms like A\ only coupling to the modulus of the order 

parameter is possible. Considering the form of Fgj the 

(2) 

elastic constant corrections due to FI^_^q can be read 
off immediately as 

Ac^ = AcO'(l-;^) 

(28) 

therefore the anomalies due to quadratic strain coupling 
are continuous consisting in a slope change below Th 
in contrast to the discontinuous jumps caused by linear 
strain coupling. The amplitudes Ap are given by domain 
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independent expressions for the two E phases which are 
now Hnear in the coupling constants g'^'- 

Eil,l): Ac0'=4,f(^) (29) 

Since oiily the (tiny) slope changes are observed in cr(T) 
m Ref. M this may well be the dominant mechanism. 
We note however that a similar coupling as in Eq. ((27)) 
exists for the nondegenerate A2, Si, B2 HO phases. Sim- 
ply + 0^) has to be replaced by (j)"^ and the slope 
change anomalies are the same as in Eqs. (|28l29p except 
for an overall factor of two. Therefore we conclude that 
the elastic constant anomalies of the type observed in 
Ref. [l3 may not be very useful for discrimination between 
the possible symmetries of HO. On the other hand these 
measurements were still done in samples with parasitic 
small moment AF phase which may obscure the results. 
Therefore they should be repeated using the pure HO 
samples. 

VI. CONCLUSION 

In this work we have reconsidered the question of 
proper symmetry of the hidden order in URu2Si2. A 
phenomenological Landau analysis based on D^h irre- 
ducible HO parameters has been applied to recent os- 
cillatory torque experiments in rotating magnetic field. 
We have concluded that the experimental observations of 
twofold torque oscillations are only compatible with the- 
oretical predictions from the two component (antiferro-) 
quadrupolar E(l,l) {Oyz,Ozx) phase . The second sin- 
gle component phase E(1,0) can be ruled out because of 



a shifted angular and wrong temperature dependence of 
the twofold amplitude. For single domain samples of the 
E(l,l) phase exhibits twofold oscillations with the proper 
field angle and temperature dependence. 

All one dimensional representations, namely 
quadrupolar Bi{Ox2-^y2) and B2{Oxy) as well as 
hexadecapolar A2{Hxy[x'^-y2-^) are incompatible with 
torque experiments. We conclude that the two com- 
ponent E-type quadrupole is the HO symmetry in 
URu2Si2. This order parameter is also compatible 
with observations in uniaxial stress experiments^^ and 
resistivity anisotropy^S. The elastic constant anomalies 
observed sofar are possible in the E type hidden order 
but may presently not be used for discrimination of HO 
symmetry. 

It remains to be seen whether the quadrupolar E(l,l) 
type HO can be identified directly. In this respect res- 
onant x-ray scattering is a powerful method. Results 
obtained before^ exclude the Bi{Ox2-y2) and B2{Oxy) 
type quadrupoles in agreement with present torque re- 
sults. However possibility of quadrupolar E(l,l) type 
HO has sofar not been investigated by this method. One 
should be aware that due to the domain problem of the 
degenerate quadrupole it could only be used for small 
single domain samples. 
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